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We propose to realize a two-dimensional chiral topological superconducting (TSC) state from the 
quantum anomalous Hall plateau transition in a magnetic topological insulator thin film through 
the proximity effect to a conventional s-wave superconductor. This state has a full pairing gap in 
the bulk and a single chiral Majorana mode at the edge. The optimal condition for realizing such 
chiral TSC is to have inequivalent superconducting pairing amplitudes on top and bottom surfaces 
of the doped magnetic topological insulator. We further propose several transport experiments to 
detect the chiral TSC. One unique signature is that the conductance will be quantized into a half¬ 
integer plateau at the coercive field in this hybrid system. In particular, with the point contact 
formed by a superconducting junction, the conductance oscillates between e^/2h and e^/h with the 
frequency determined by the voltage across the junction. We close by discussing the feasibility of 
these experimental proposals. 

PACS numbers: 74.45.-l-c 73.43.-f 71.10.Pm 73.40.-c 


I. INTRODUCTION 

The search for topological states of matter has become 
a central focus in condensed matter physics. Chiral topo¬ 
logical superconductors (TSC) in two-dimensions (2D) 
with an odd-integer Chern number are predicted to host 
a Majorana zero mo de i n the vortex core, which obeys 
non-Abelian statisticJ^^ and has potential applications 
in topological quantum computatioiP. A chiral TSC with 
Chern number Af breaks time-reversal symmetry, and has 
a full pairing bulk gap and N' topologically protected gap¬ 
less chiral Majorana edge modes (CMEMs), which can 
be viewed as a superconducting analogy of the quantum 
Hall (QH) state®£l. As a minimal topological state in 
2D, the N = 1 chiral TSC is of particular interest, as 
its edge state has only half the degrees of freedom of the 
QH state with Chern number C = 1. Intensive efforts 
have been made to search for the chiral TSC in 2EPHIZ1, 
however, it has not yet been confirmed in experiments. 

In principle, a QH state with Chern number C in prox¬ 
imity with an s-wave superconductor (SC) can be natu¬ 
rally viewed as a chiral TSC with even number Af = 2C 
CMEMs. Therefore, it is theoretically possible to re¬ 
alize a chiral TSC with odd number of CMEMs near 
a QH plateau transitioiP^. However, the strong mag¬ 
netic field required in a QH state will severely hinder 
the superconducting proximity. Instead, the quantum 
anomalous Hall (QAH) state has a finite C hern n umber 
C in the absence of an external magnetic fielcP^E^, which 
has been theoretically predicted in magnetic topological 
insulators (TIs) with ferromagnetic (EM) ordering^SlEni 
and experimentally realized (for C = ±1) in both Cr- 
dopecPimil and V-dopecP^ (Bi,Sb) 2 Te 3 magnetic TI thin 
films. More recently, a new zero-plateau QAH state with 
C = 0 and the plateau transitions among C — ±1,0 states 
have been theoretically predictecP^ and experimentally 


observecP^^^. Without requiring a large external mag¬ 
netic field, the plateau transition from the C = ±1 QAH 
to the zero-plateau C = 0 state is a unique parent system 
for realizing a Af = ±1 chiral TSC. 

In this paper, we propose to realize the A/” = ±1 chiral 
TSC in a magnetic TI near the QAH plateau transition 
via the proximity effect to an s-wave SC. The optimal 
condition for realizing the chiral TSC is to have inequiv¬ 
alent SC pairing amplitudes on top and bottom surfaces 
of the doped magnetic TI. We then propose several trans¬ 
port experiments to detect this chiral TSC. Generally, 
the conductance could be quantized into a half-integer 
plateau at the coercive field in this hybrid system (Fig.[^, 
as a signature of the neutral CMEM backscattering. In 
particular, with a point contact formed by a SC junction 
(Fig. 1^, the conductance oscillates with a frequency de¬ 
termined by the voltage across the junction. Lastly, we 
briefly discuss the temperature dependence on the trans¬ 
mission of CMEM and the feasibility of these experimen¬ 
tal proposals. 

The organization of this paper is as follows. After this 
introductory section, Sec. H describes the effective model 
for the SC proximity effect of the QAH state in a mag¬ 
netic TI thin film. Section HI presents the results on the 
phase diagram, edge transport and experimental propos¬ 
als on point contacts. Section IV presents discussion on 
the feasibility of experimental realization of chiral TSC 
in a magnetic TI. Section V concludes this paper. Some 
auxiliary materials are relegated to appendixes. 


II. MODEL 

To start, we consider the SC proximity effect of the 
QAH state in a magnetic TI thin film with EM order. 
Without the proximity effect, the low energy physics 
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FIG. 1. (color online). The hybrid QAH-SC device. In region 
II, a chiral TSC state is induced through the proximity effect 
to an s-wave SC layer on top of the QAH in magnetic TI. A 
back-gate voltage Vbg is applied to control the Fermi level in 
region II. Voltages Vi and V 2 are applied on leads 1 and 2, 
respectively. The SC layer is grounded through a lead in its 
bulk. 

of the system only consists of the Dirac-type surface 
states (SS)P^. The 2D effective Hamiltonian is Ho — 
Ek^k^o(k)V'k, with V'k = and 

iJo(k) = kyaj;T:, - kxCTyTz + rn{k)T^ + Xa^, (1) 


where H±{k) = kya^ T k^dyC^z + {m{k) ± X)az<,z T AcTyCj/ 
with <;x,y,z the Pauli matrices in Nambu space. The topo¬ 
logical property of is clearly seen by a further basis 
transformation into a block diagonal form: 

, (4) 

where h±{'k) = kyUx — k^Uy + {m{k) -I- A ± |A|)tT 2 char¬ 
acterizes a Px H ipy SCPJ^. The BdG Chern number of 
h± (k) depends only on the sign of mass m{k) -I- A± | A| at 
the r poinlP^. Therefore, the Chern number of 7J+(k) is 
A/+ = —2 for |A| < — toq — A, A/+ = —1 for |A| > |too + A| 
and A/+ = 0 for |A| < toq + A. Similarly, the Chern num¬ 
ber of if-(k) is N- = 2 for |A| < A — mo, N- = 1 for 
|A| > |mo — A| and N- = 0 for |A| < mo — A. The to¬ 
tal Chern number of the system is then N = A/+ -I- N -. 
Fig. 1^ shows the phase diagram of the system. The 
phase boundaries are determined by A ± (mo ± A) = 0, 
which reduce to the critical points A = ±|mo| between 
the C = ±1 QAH and the zero plateau normal insula¬ 
tor (NI) for A = 0. An infinitesimal SC gap drives the 
QAH phase into a N = ±2 TSC. More importantly, the 
A/” = ±1 TSC state emerges in the neighborhood of the 
transition between the QAH phase and NI phase. 


where Cko- annihilates an electron of momentum k and 
spin cr and superscripts t and b denote SS in the 

top and bottom layers, respectively, at and t) {i = x, y, z) 
are Pauli matrices for spin and layer, respectively. A is 
the exchange field along 2 : axis induced by the FM or¬ 
dering. Here A oc (S') with (S) being the mean field 
expectation value of the local spin, and the value of A 
can be changed during the magnetization reversal pro¬ 
cess in magnetic TIs. m{k) = mo + mi(A;^ + ky) de¬ 
scribes the hybridization between the top and bottom 
SS. The Chern number of the system is C = A/|A| for 
|A| > |mo|, and C = 0 for |A| < |mo|. Correspond¬ 
ingly, the system has \C\ chiral edge stated. In proximity 
to an s-wave SC, a finite pairing amplitude is induced 
in the QAH system. The Bogoliubov-de Cennes (BdC) 
Hamiltonian becomes ’Hbug = X^k'^^k^Bdc'I'k/S, where 
= [(4^, 4;. 4 t 4;). (c-kt: 44;, 44;, 44;)]^ and 


H^dG = 

Ak = 


fHo{k)-pi Ak 

A], -iJ*(-k) + /r 

/ ilXlCFy 0 \ 

0 il^2dy) ' 


( 2 ) 


Here p, is chemical potential, Ai and A 2 are pairing gap 
functions on top and bottom SS, respectively. 

In a simple case for /r = 0 and Ai = —A 2 = A, a 
basis transformatiorP^ decouples the BdC Hamiltonian 
into two models with opposite chirality, and 


(H+{k) 0 A 

V 0 H^{k))' 


( 3 ) 


III. RESULTS 
A. Phase diagram 

Now we turn to the optimal condition for realizing the 
N = ±1 TSC. First, consider the phase diagram for 
p = 0 and general values of Ai and A 2 . The phase 
boundaries are determined by the bulk BdC gap closing 
in Eq. Q. Assuming A 2 = oAi and a is real, the phase 
boundaries are given by =F(l — a)AiA-|-A^ = rrig-l-aA^, as 
shown in Fig. For Ai = A 2 , the Chern number jumps 


Ai=A, A2=0, |J=0 


(b) 



A,=A2=A, |J=0 



FIG. 2. Phase diagram of the QAH-SC hybrid system with 
typical parameters, (a) Ai = A, A 2 = 0, p = 0. (b) Ai = 
A 2 = A, p = 0. (c) Ai = A, Aa = 0, /r = 0.7. (d) Ai = 
— A 2 = A, p = 0.7. Here Ai, A 2 , p are in the units of |mo|. 


^^BdG = 
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directly from M = ±2 to A/" = 0, and A/” = ±1 TSC 
phases disappear due to accidental particle-hole symme¬ 
try in Hq with fi = 0. As A 2 decreases, the Af = ±1 
TSC phase space emerges and becomes the widest at 
A 2 = 0. In particular when A 1 A 2 < 0, a helical TSC 
phase with helical Majorana edge states emerges on the 
A = 0 line (Fig. 1^). The general case for complex a 
is studied in Appendix A, where the topology of phase 
diagram remains unchanged. Next, for the case ^ 7 ^ 0, 
which corresponds to the SC proximity effect of a doped 
or electrically gated QAH system, the proximity effect is 
effectively enhanced by the finite density of states at the 
Fermi leveP^. As shown in Fig. the phase space of 
A/” = ±1 TSC near the A = 0 axis enlarges from p, = 0 
to /i 7 ^ 0. Therefore, the optimal condition for Af = ±1 
TSC is /r 7 ^ 0 and A 2 = 0. This leads us to design 
the transport device in Fig. The s-wave SC is only 
grown on top of the magnetic TI in region II to ensure 
the proximity pairing gap of the top SS is larger than that 
of the bottom SS, while the Fermi level can be tuned by 
the back-gate. The size of the SC layer should be larger 
than the back-gate electrode so that there is no metal¬ 
lic regions in the device. Similarly, one can also employ 
another device geometry by using a global back-gate and 
two top-gates in region I and III, to tune the Fermi levels 
in region I, II, and III separately. 


B. Edge transport and half-plateau 

To identify the A/” = 1 TSC in the QAH-SC hybrid 
system, one can probe the neutral Majorana nature of 
CMEM or trap the vortex core zero mode. Several 
methods have been proposed to measure the Majorana 
fermionsPSHm. Here, we base our discussion on a recent 
proposal studying the CMEM backscatterinj^. The ba¬ 
sic setup is shown in Fig.[l] consisting of a magnetic TI in 
proximity with a grounded top SC layer in region II and 
two current leads at the corners. When the magnetic do¬ 
mains of magnetic TI are aligned in the same direction, 
the magnetic TI is in a QAH state with a single chi¬ 
ral edge state propagating along the sample boundary. 
During the flipping of the magnetic domains at the co¬ 
ercive field, A decreases and the magnetic TI enters the 
NI with a zero-plateau in Hall conductance Uxy over a 
finite range of magnetic fielcpZlMl^ g^g shown in Fig. §3. 
Either perpendicular or in plane external magnetic field 
could induce such plateau transitiorP^. When the SC 
proximity effect is sufficiently strong, the superconduct¬ 
ing region H experiences the BdG Chern number varia¬ 
tion A/’=—2—>•—1—>-0—>-1—?'2asA decreases in the 
hysteresis loop (dashed line in Fig. [^). Therefore, the 
transport setup Fig.j^is a QAH/NI-TSC/NSC-QAH/NI 
junction. As we will discuss in details below, the edge 
transport features of the junction uniquely convey the 
topological properties of the SC in region H. 

The QAH edge state can be viewed as two CMEMs 
since a C = 1 QAH state is topologically equivalent 


to a A/” = 2 TSC. Therefore, in the case of QAHc=i- 
TSCa/'= 2 -QAHc=i junction (Fig. |^), the edge current 
will be perfectly transmitted. By contrast, if the junction 
is QAHc=i-TSC 7 V=i-QAHc=i (Fig. [^), the chiral edge 
state in the QAH regio n separates into two CMEMs at 
the TSC boundarjP^. One CMEM is perfectly trans¬ 
mitted, while the other is totally reflected. The edge 
transport of the junction is governed by the generalized 
Landauer-Biittiker formalism, which includes the contri¬ 
butions f rom b oth the normal scattering and Andreev 
scatterinj^^m^ Hhe general relationship between current 
and voltage on lead 1 and 2 shown in Fig. is /i = 
(eVh)[(I - 7^ + 7 ^ 2 l)(Vi - 140 ) - {V - rxm - K?)], and 
h = (eV/l)[(l-7^' + 7^4)(F2-F,?)-(r-Xl)(Fl-I4?)]. 
Here 14° = 0 is the voltage of the grounded SC layer, Ji 
and I 2 are currents flowing into leads 1 and 2 , respec¬ 
tively. TZ, T, TZa and Ta are the normal reflection, nor¬ 
mal transmission, Andreev reflection and Andreev trans- 



FIG. 3. (color online), (a) Phase diagram of the QAH- 
SC hybrid system for g = 0 and Ai = —A 2 = A. Only 
A > 0 is shown, (b) Without SC proximity effect, the 
Oxy = — 1 —>■ 0 —>■ 1 QAH plateau transition occurs at the 
coercivity when the magnetization flips, (c) With SC prox¬ 
imity effect to region H in hybrid device Fig. (712 shows 
plateau transition 1—>1/2—>0—>1/2—^-lin the hysteresis 
loop. The half-integer plateau in (J 12 manifests the N = 1 
TSC. (d)-(j) The edge transport configuration at A, B, C, D, 
C/ B' and A' in (c). There is no backscattering for Af = ±2 
TSC in (d),(j), and Majorana backscattering for A/” = ±1 TSC 
in (e),(i). Red and blue arrows represent (c ± cQ CMEMs, 
respectively. NSC: normal, topologically trivial SC. 
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mission probabilities for an electron injected from the 
left, while TZ', T', and 7^ are for an electron coming 
from the right. The two-terminal conductance is then 
defined as ai 2 = I/{Vi - V 2 ) = (Ii - /2)/2(Vi - ¥ 2 ). 
For the QAHc=i-TSCa/'=i-QAHc=i junction in Fig. [^, 
the probabilities of normal scattering and Andreev scat¬ 
tering are equaP^, and we have TZ = TZa = T = Ta = 
TZ' = TZ'j^ =T' = Ta = ^/‘^, resulting in a half-quantized 
conductance 

Moreover, since the SC layer is not floating but grounded, 
the quantized net incoming current /sc = (Vi + V2)e'^(h 
will be flowing from the SC layer to ground. Here we 
point out that the supercurrent due to the phase fluctu¬ 
ation of SC order parameter may give a small correction 
to conductance, which scales as where i is the 

width of CMEM, and L is the size of SC. For an esti¬ 
mation, / ~ 0.5 /im, therefore such correction is neglible 
for L > 50 /im. In contrast, the N = 2 TSC junction in 
Fig. exhibits a quantized conductance (7\2 ~ eV/P. 

The entire plateau transition of <712 in the hybrid junc¬ 
tion device is shown in Fig. In correspondence to 
the QAH plateau transition of axy in Fig. (Ti 2 also 
exhibits plateaus quantized at jh and 0 when region 
II is Af = ±2 TSC and Af = 0 NSC, respectively. 
In addition, an intermediate half-quantized plateau at 
e^/2/i could occur at the coercivity under the condition 
|A|-|-|mo| > |A| > I mo I, which is a unique signature of the 
A/” = ±1 TSC in region II. We emphasize that a plateau 
usually indicates a stable phase instead of a fine-tuned 
state. The size of backscattering region is not necessar¬ 
ily mesoscopic. In fact, the size L of the TSC region 
sets a temperature scale ksTint ~ vm/L, above which 
the interference effect vanishes due to thermal averaging, 
where vm is the Fermi velocity of CMEM. For an esti¬ 
mation, L ~ 200 /tm, Vm ~ 2.0 eV A, ~ 10 mK. 
Therefore, the half-plateau is robust at large L and finite 
temperature T > Tint- The plateau transitions and cor¬ 
responding edge transport configuration in the hysteresis 
loop are illustrated in Fig. |^-j. In particular, four 1/2- 
plateaus occur around the critical magnetic fields zLH* 
and ±//| shown in Fig. [^. 


C. Point contact 

Another useful transport configuration is a point con¬ 
tact formed by two SC islands which allow the transmis¬ 
sion of CMEMs, as shown in Fig. |^. A voltage I4c is 
applied onto island TSCi, while TSC 2 is grounded. If 
either TSCi or TSC 2 is a Af = 2 TSC, the edge current 
will be perfectly transmitted. Non-trivial physics occurs 
when both TSCi and TSC 2 are Af = 1 TSC. An incident 
edge electron from bi splits into two CMEMs, one is per¬ 
fectly transmitted along the edge, while the other is scat¬ 
tered at the point contact with transmission amplitude t. 



FIG. 4. (a) The point contact configuration of two SC is¬ 
lands with SC phases 4>i and (f> 2 , across which the reflection 
and transmission amplitudes of the CMEMs are r and t. (b) 
The conductance cti 2 as a function of r for different coupling 
strengths C A dc current Hows between ai and 02 , an ac volt¬ 
age between them is measured, with frequency / = 2eVsc/h. 


which depends on the phase difference 5(j) = of two 

TSCs (see Appendix B). The I-V relation in this geome¬ 
try is/i = {eyh)[{l-TZ + TZA){Vi-V,n)-{r-U)V 2 ], 
and I 2 = (eVh)[(l -TZ' + TZ'jf)V 2 - {T-Ta){Vi - Esc)], 
where TZ = TZa = TZ' = TZ'a = rV4, T = V = {l + tf/A, 
Ta = Ta = (T — t)^/4, r is reflection amplitude and 
-I- = 1. Therefore, I = e^{l + t){Vi — V 2 — 14c)/2h. 

Note that the current is proportional to the tunneling 
amplitude t, not the tunneling probability. If Ec = 0, we 
have CT 12 = {l + t)e^/2h, which directly measures t of the 
neutral CMEMs. A finite 14c leads to a time dependent 
(5(/, which in turn affects t. A simple tunneling model for 
the CMEM is (also see Appendix B) 

^ftunnel = iCTzOx - k{x) Sm{S<p/2 - 4>o)<Xy, (6) 

where k{x) is nonzero in a finite interval, and the basis 
is the CMEMs ( 71 , 72 ) shown in Fig. The transmis¬ 
sion amplitude t at zero-energy in this model is t((5(/) = 
1/cosh[^sin((5(()/2 — (/o)], where C = / dxK{x)/2. Within 
this model, t is purely real. With a fixed Esc across the 
point contact, 54> varies linearly with time t with a slope 
dScjj/dT = 2eVsc/k. We can define a new conductance 

which is a periodic function in time with the Josephson 
junction frequency / = 2eEc/h. Fig. & shows 17(2 a 
function of time for different values of f. The time oscilla¬ 
tion shape of a '12 are different for a weakly coupled point 
contact (small f) and a strongly coupled one (large ^). 
However, a '12 always oscillates between e^/2/i and jh, 
since there is always at least one perfectly transmitted 
CMEM, which is also a unique feature of the A/” = 1 
TSC state. 


D. Temperature dependence 

We further consider the temperature dependence of 
the above CMEM transmission (see Appendix D). It is 
straightforward to see by a dimensional counting that 
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in the above free Majorana fermion model is 
marginal, therefore it remains constant at low temper¬ 
ature T. When the leading four-fermion interaction (ir¬ 
relevant) is included, the tunneling amplitude acquires 
a weak temperature dependence. For 14c = 0, in this 
case 17(2 = <^121 the renormalization group analysis gives 
a power-law correction St ~ to t, where Xp is 

the bare fermion interaction strength. The conductance 
a '12 oc (1 + t) will therefore decrease as T increases. This 
perturbative result is no longer valid above a characteris- 

_ 1 /Q 

tic temperature of Tc ~ Ap ' , when the correction St is 
comparable to t. For higher temperature Tc < T \A\, 
t will flow towards 0, and the two TSC islands will behave 
like a single connected TSC analogous to that shown in 
Fig.g In this regime, one can formulate a similar point- 
contact tunneling model between the left and right edges 
of the new TSC as in Eq. ([^, but with an additional 
vortex tunneling through the bulk TSC. At high tem¬ 
perature, the leading contribution to t then comes from 
the vortex tunneling, which leads to t ~ where 

Ao- is the bare vortex tunneling strength. Therefore, the 
half-quantized plateau in (T 12 remains robust in the high 
temperature regime Tc < T <C |A|. 


could be quantized into a half-integer plateau at the co¬ 
ercive field in this hybrid system, as a unique signature of 
the neutral CMEM backscattering. We emphasize that 
such an experiment can work at reasonable temperature 
and does not depend on the interference effect of CMEM. 
We hope the theoretical work here can aid the search for 
chiral TSC phases in hybrid systems. 
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Appendix A: Phase diagram under complex 

a = A2 / Ai 


IV. DISCUSSION AND EXPERIMENTAL 
REALIZATION 

Einally, we discuss the feasibility of our proposals. 
Experimentally, to observe the A/” = ±1 chiral TSC 
and all of the four half-quantized conductance plateaus, 
a good proximity effect between SC and magnetic TI 
is necessary. Moreover, the critical field of SC 
should be larger than the coercivity Hi 2 in magnetic 
TI. From Ref. [3H| and |3S1 the estimated H* ^ 0.05 T 
and H 2 ~ 0.2 T. The candidate SC materials are Nb and 
NbSe 2 . The bulk Nb is a type I SC with Tgc = 9.6 K and 
^ 0.2 T, while a thin film Nb becomes a type II SC 
with upper critical field ~ 1 T. NbSe 2 is a type II 
SC and shows good proximity effect with Bi 2 SeaP^even 
at 4.2 K and 0.4 T, where the proximity effect induced 
SC gap is A ^ 0.5 meV. The width of the CMEM i can 
be estimated as vp/A ~ 0.52 /rm, where the Fermi ve¬ 
locity vp ~ 2.6 eV For a typical junction voltage 
I4c 1 fXV, f ^ 0.48 GHz, which is easily accessible in 
experiments. 


V. CONCLUSION 

In summary, we propose to realize the A/” = ±I chiral 
TSC in a magnetic TI near the QAH plateau transition 
via the proximity effect to an s-wave SC. We show that 
inequivalent SC pairing amplitude on top and bottom 
surfaces in doped magnetic TIs will optimize the A/” = ±I 
chiral TSC phases. Several edge transport measurements 
have been proposed to identify such J\f = 1 TSC in the 
QAH-SC hybrid system. In particular, the conductance 


In the paper we have only considered the case a = 
A 2 /A 1 is real. In general, in the absence of time re¬ 
versal symmetry (as is in our model), a = |a|e*'^“ is 
complex. Correspondingly, the phase diagram will be 
modified quantitatively, but the topology of the phase 
boundaries remains unchanged compared to those shown 
in Fig. of the paper. 

By a proper choice of basis we can always set 
Ai = A real. As an illustrative example, we con¬ 
sider here the case |a| = 1, namely a = A 2 /A 1 = 
Via a unitary transformation —>■ 

^2 is transformed into 
a real number A 2 = Ai = A, while the hy¬ 
bridization m{k) between the top and bottom SS be¬ 
comes a complex number e“*'^“'^^m(/c). Therefore, we 
can always set two of the three parameters Ai, A 2 
and Too to real numbers. Diagonalizing the BdG 
Hamiltonian iFedC yields the energy spectrum = 


X±H[m{k) sin((()„/2) ± A] -|- mikY cos^ (Y>^j2) 


The phase boundaries are given by the gap closing of the 
energy spectrum: 


A ± (too sin((^Q,/2) ± A)2 -|- tuq cos^ ((^q,/ 2) = 0, (AI) 
namely, the following hyperbolas: 

A^ — (a ± mo sin((()o,/2)^ = m^ cos^ (<(>a/2). (A2) 

The phase diagram is shown in Fig.[^ As one can see, the 
topology of the phase diagram does not change much. In 
particular, when (()„ = 0 and tt, the phase diagram is as 
indicated in Fig. 2b and Fig. 3a of the paper, respectively. 
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FIG. 5. The phase diagram for A 2 = e*'^“Ai and ^ = 0. 
When 4>a = 0, the M = ±1 TSC phases disappear, while 
when (j)a = TT, the phase spaces of A/” = 1 and Af = —1 TSC 
touch each other, as indicated in Fig. 2b and Fig. 3a of the 
main text, respectively. 


Appendix B: Derivation of the effective tunneling 
Hamiltonian 


Without loss of generality, consider the case |A| > A — 
mo > 0. The QAH has Chern number C = 1 and the 
SC in region II has BdG Chern number JV — 1, both 
of which come from the lower block iJ_(k) of the BdG 
Hamiltonian i?BdG- When the pairing amplitude of the 
superconductor is A = |A|e*‘^ = Ai = — A 2 with a phase 
(k) can be rewritten as 


iJ_(k) = 


h'Jk) = 


fh'^k) 0 \ 

0 -h'l{-k))’ 

/m(fc) —A±|A| —ikx^ky \ 
y ikx ± ky —m(A:) + A =F I A|y ’ 


(Bl) 

(B2) 


under the following new basis -^(e 


used the notation 




Ck^ 


Ckt + ), where we have 


--k^’ 


Ck^ 


Ckt = 


r* — 

Ckt ^kt 


V2 


and 


cii + cf 


Ckt = 


kt 


V2 


(B3) 


(B4) 


The Majorana edge state between the QAH (where | A| = 
0) and the TSC (where |A| > A — mo > 0) is given by 
h'^{k). 

As shown in Fig. [^of the paper, the lower TSCi and 
the upper TSC 2 have superconducting phases and (/)2 
respectively. For simplicity, we shall approximate m(k) 
as mo, which does not change the topological physics. If 
the upper edge of the lower TSCi is set as y = 0, the 
Hamiltonian of the corresponding Majorana edge state 
can be derived as 


Hi 


= y da; i'yi{x)dx'^i{x), 


(B5) 


where 


71 (a;) = 


g-z0i/2^i(a.) ^gi0i/2^t(a.) 


Cl (a;) = J 


e 

—00 

— 27r/4 


V2 

(\A\e>{-y)+mo-\)y 


t/4 




+e 


ci{x,y) dy, 


(B 6 ) 


(B7) 


with 0(y) defined as the Heaviside function. Similarly, 
the lower edge of the upper TSC 2 at y = yo > 0 has a 
low energy Hamiltonian 


where 

72 ( 3 ;) = 
C 2 (a;) = 


H 2 = - j dx ij 2 {x)dx^ 2 {x), 


(B 8 ) 


72 

„(A-mo-|A|e(i/-i/o))!/ 


(B9) 


^—i'K/4 


ct(a;,y) 


+e“/7;(x,y) 


dy. 


(BIO) 


We shall assume the point contact extends in the interval 
0 < a; < L, and the two edges have a nonzero hopping 
and pairing term: 

i-L 


Hi = - 


dx 


Jhc\ {x)c 2 ix) 


+Jp(Aj; + A 2 )ci{x)c 2 {x) + h.c.], (Bll) 

where A 1.2 = |A|e*‘^c 2 ^ When projected into the low 
energy Hilbert space of 71 and 72 via the substitutions 


Cle®'^^/^7i/72, C 2 -)► e 


* 02/2 


72 


/V2, 


(B12) 


this term becomes: 

I'L 

Hi = 2 


J dx iK(x) sin - 4>o^ 71 ( 3 ^) 72 ( 2 ;) 


= 2 


dx iX{x)ji{x)^ 2 {x), 


where 


6(j) = 4>i- 4)2, 
k { x ) = \ Jh/‘2. + Hm(Jp)|, 

4)0 = arg [Jh + i 2 Im( Jp)] . 

For simplicity we have defined 

A(a;) = k { x ) sin f y “ ' 


(B13) 


(B14) 

(B15) 

(B16) 


(B17) 


The total tunneling Hamiltonian is then Tftunnei = Hi + 
H 2 + Hi as given in Eq. (§ of the paper. The eigenwave- 
function 1 // = (yi, 772 )^ at energy E can then be obtained 
by solving the following Shrodinger equation: 


/5a; /A(a;)A /ryi 

-i\{x) -idx J \r ]2 


= E 


(B18) 
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The solution for a wave incident from x = —oo with Upon basis transformation from Majorana fermions to 
momentum k is E = k, and charged fermions on QAH edges, we have 


irii{x),V2{x)) 


^—ikx A sinh(\/A^ —fc^L) 

Q[L] _ ® ; 

^ Q[L — x\ A sinh[\/ X'^ —k‘^{L — x)\ 

G[L] ’ WA 

VA 2 -fc 2 e-’'=" nA 
g[L] ’'JJ 


{x < 0 ) 

(0 < a; < L) 

{x > L) 

(B19) 


where function G[x] = -s/A^ — cosh(-\/— kP'x) — 
ik smh{y/— k'^x). At low energies k the wave- 

function can be approximately written as 


{vi{x), mix)) 


X{x')dx' 


, sinh 



Xix')dx' 


(B20) 


from which the transmission and reflection amplitudes 
can be extracted out as 


/ ai,k \ 
Oq-k 


/ r 

r t -I- 1 t — 1\ 

^Uk 


r 

1 —1 

-f 

1 

02,k 

2 

t* + 1 

t* - 1 -r* -r* 

d2,k 

V4,-k/ 


\t* - 1 

r + 1 -r* -r* / 

\bi-J 


(C2) 


based on which the normal/Andreev transmis¬ 
sion/reflection probabilities are given as T = |t + lp/4, 
Ta = |A — lp/4, and TZ = TZa = kP/4- According to the 
generalized Landauer-Biittiker formula, the conductance 
defined in the main text is 


(T 


I _ 

12 — 


1 + Re(t) 

2 ~h' 


(C3) 


Note that the conductance a [2 merely depends on the 
real part of Majorana transmission coefhcient t, phys¬ 
ically it is due to the fact that charged fermions are 
treated as combinations of Majorana fermions with trans¬ 
missions t and perfect transmission 1 . 


cosh i f dxX(x)) cosh [^sin((5(/)/2 — (/o)] ’ ^ ^ 

r = tanh dxX(x)^ = tanh sm(S(/)/2 — c/o)], 

(B22) 

where ^ = / dxnix). Note that t is always real and pos¬ 
itive at low energies. For scattering at a finite energy 
E = k, the transmission amplitude t is generally com¬ 
plex. 


Appendix C: S-matrix and conductance in general 
Josephson junction setup 

Here we formulate the scattering matrix of edge states 
in the setup of Fig. 4a, and derive the conductance a'i 2 - 
The edge fermions at four ends of the sample are denoted 
by ai .2 and 61.2 as shown in Fig. 4a. With transmission 
coefficient t and reflection coefficient r at the point con¬ 
tact, the scattering matrix S due to the point contact 
is 


Appendix D: Temperature dependence and 
renormalization group analysis 

In this section we analyze the temperature dependence 
of Majorana transmission coeff icient t by renormalization 
group technique in detaiP^*^. Specifically, we focus on 
its real part Re(t), since it is directly related to the con¬ 
ductance a 12 - Our starting point is the action for the 
model in Eq. of the paper, 

5o = y di- y dx[-fiiidr + 9x)7i + 72*(i9r - dx)l2 
-I-2^<5(x)sin(i5(//2 - (/o)i7i72]- (Dl) 

Since the Majorana tunneling occurs locally at x = 0, 
the scaling dimension of the tunneling strength ^ van¬ 
ishes, i.e. [^] = 0. Therefore, ^ is invariant when the 
temperature T of the system changes, and so does the 
transmission coefficient t. 

The temperature dependence of t comes from higher 
irrelevant terms at the point contact. The leading irrel¬ 
evant term is a four fermion interaction of the following 
form: 


/ fli.k + o, 
fli.k — a 
02 ,k + a 


= S 


fbiM + ^i,-k\ 

k — 

^ 2 ,k + bt _] 


2,k 

®L-k/ 


V^2,k ^2 _i^ 

/r 

0 

t 


'^^'l,k + ^l,-k\ 

^l,k ~ ^1,-k 

0 

0 

0 

1 1 

t* 

0 -r* 

0 

^'2,k + bl 

Vo 

1 

0 

0/ 

\b2.M-bi-J 


Hp = J dxXp5{x)^idxmi2dxl2- (D 2 ) 

It represents the tunneling of one pair of Majorana 
fermions from one edge to the other. The scaling dimen¬ 
sion of Xp is [Ap] = —3, hence it is irrelevant and scales as 
Ap® ~ XpT^ when T ^ Q. Increasing the temperature T 
will enhance the effective interaction strength Ap, which 
affects the transmission coefficient t. 

The contribution of Edp to the transmission coefficient 
t can be calculated perturbatively as follows. Suppose 

































both ^ and Xp are small, so that perturbation theory can 
be used. We shall regard Hj = 2i^S{x) sm{6(j)/2—(j)o)jij2 
and Hp given above as the perturbation. Consider an 
in-state |i) = of Majorana fermion 71 , and a 

transmitted out-state |/) = 7 i,_/c'|n), where |n) is the 
system ground state. The transmission coefficient t is 
then given by 

t « (/I |i) , (D3) 

where stands for the time ordering. The zero-order 
is simply 5kk'■ The first-order contribution is 

= (Hi + Hp)dr |i) . (D4) 

Since Hj is odd in 71 and 72 , its first-order contribution 
vanishes. The second term of Hp is purely imaginary 
and therefore does not contribute to the conductance ( 7 ^ 2 ■ 
The second-order correction 

~ (/I Tr I {Hi + Hp){t){Hi + Hp){T')dTdT' |z) . 

(D5) 

The H] term gives a constant contribution ^ 
—sin^(5(/)/2 — ^o)j in agreement with calculations in 
Appendix B. The cross term HjHp vanishes because it 
is odd in 71 and 72 . The H^ term results in a tempera¬ 
ture dependent correction to the real part of transmission 
coefficient t as 

^Re(t) ^ -Skk' (Af )" = -Skk'XlT^ (D 6 ) 

Therefore, the transmission coefficient t generically de¬ 
creases as temperature T increases. When the tempera- 
ture T is above a characteristic temperature Tc ^ Xp ' , 
the interaction Xp at the point contact dominates, so that 
t becomes small and r becomes large. In this case, the 


above perturbative treatment is no longer valid. How¬ 
ever, this case can be effectively viewed as a breaking up 
of original Majorana edge states 71 and 72 and a remerge 
of them into two new Majorana edge states "01 and tjj 2 
on the left and right of the point contact, and of the 
two TSCs merging into a single TSC. In the temperature 
range <C T ^ |A|, we can do a perturbation calcula¬ 
tion about the high temperature fixed point before the 
superconducting phase is destroyed. 

This scenario is very similar with our setup in Fig. la, 
except that the two edges are brought together at the 
point contact. Since the region between the edges in this 
case is a SC, there are both fermion tunnelings and vortex 
tunnelings between edge^^. The effective action for this 
point contact is 

S' = J dr J dy [ipiiidr + Vmdy)^^! + 4’2i{dT - Vmdy)i ’2 

+X^S{y)iij}iijj2 + XaS{y)aia2 ], (D7) 

where ci and cr 2 are the vortex operators on edges with a 
scaling dimension [ui] = [a 2 ] = 1/16. Dimension count¬ 
ing renders [A,/,] = 0 and [Ao-] = 7/8, so the vortex-vortex 
tunneling is the most relevant. Therefore, at a high tem¬ 
perature T, the vortex-vortex tunneling term gives the 
temperature dependence of transmission coefficient t 

t - A^T-^/^. (D 8 ) 

The power-law relation is valid above a characteristic 
temperature T/ ^ A®^^, provided the SC gap |A| is much 
higher. In fact, this confirms the robustness of the half- 
quantized plateau. For in the setup with reasonable finite 
temperature, the edges are far away from each other, so 
the tunneling strengths including X^ are sufficiently tiny, 
resulting in an extremely low T/. 
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